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Abstract 

The periods after the end of the "primary" nucleation are consid- 
ered. The approximate analytical description is given. The process is 
split into several periods which form the loop of evolution. 



1 Pleminary remarks 

The power of metastability in the system is described by the value of super- 
saturation ( defined as the ratio of surplus number density of molecules to 
the saturated number density of molecules. After the end of nucleation one 
can show that the size droplet spectrum is rather well localized in the space 
of sizes. If we choose the free molecular regime of growth (in order to escape 
from the problems of the density profiles around droplets) it is natural to 
come to the natural size coordinate 

p = v 1 ' 3 

where v is the number of molecules inside the droplet. The law of regular 
growth is rather simple 

dt t p 

with the characteristic time T clS cL parameter and the critical size p c which 
is determined as 

p c = 2a/31n(C + l) 
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where a is the renormalized surface tension. For the small supersaturations 

p c fa 2a/3( 

On the base of the critical size the value of supersaturation can be expressed 
as 

C = 2a/3p c 

The asymptotic rate of growth for the supercritical droplets is very simple 

A C 
[ dt )as r 

From the theory of nucleation it is known that the distribution function 
is well localized in the scale of sizes. Namely, the half-width Aix in p-axis 
satisfies the following estimate 

Here /* is the amplitude of spectrum, is the "quantity" of surplus vapor 
calculated in terms of supersaturation, parameter T can be estimated as 



r ~ v, 



c 



and the value v c is the number of molecules in the critical embryo. 
The number of droplets can be estimated as 

N = f*A lX 

At the final of the surplus vapor exhaustion one can write 

Nz% n = *. 

where Zfi n is an imaginary final coordinate in p-scale (actually there is no 
real final coordinate it is some quasistationary approximation). 

We consider the situation of decay here but the situation of dynamic 
conditions can be easily considered by quasistationary approximations and 
by the first step of corresponding iteration procedures. 

Then we come to 

^HL = r 1 / 3 ~ z// 3 (1) 
Atf 

Since y c > 1 we see the certain hierarchy. 
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2 Regular relaxation of the size spectrum 



When the asymptotic law of the regular growth is used the evolution of the 
system after the end of nucleation [I] can be described by means of the first 
order differential equation 



with "coefficients" as the full momentums of the distribution function 



Here f(x) is the distribution density determined by the nucleation theory [I]. 
The variable a; is the shift at p-scale of the size of droplet from the "marked" 
size z of a droplet born at some moment. 

If we choose as z the size of the droplet formed at the initial moment of 
time then _ 



The initial condition for the differential equation is the following: At 
t — t* + Ait the position of the size spectrum is z = A\x. Here t* is the 
"characteristic" moment, Ait is the duration of the "back" side spectrum 
formation. Approximately Ait = Aixr/§*. So, t* + Ait is the moment 
of the "end" of the nucleation period. Also here Aix is the length of the 
back side of the size spectrum. All these characteristics can be given by the 
nucleation theory. 

This differential equation can be easily integrated 



dz 
~dt 






dz 



Condition of the end of the regular relaxation is 



(2 - 3)p c = z-Aix 



(2) 



or 



(2-3)p, 



'c 



= Z 



With appropriate renormalization. 
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Since 

2a 2a 
Pc = 31n(C + l) ~ 3C 

and 

i 

it is an ordinary algebraic equation which can be effectively solved on the 
base of the strong inequality 

At the moment when the equation (0) is satisfied one has to take into 
account the diffusion correction to the law of droplets growth. 

The final values can be made more accurate. We can use the balance 
equation 

$ = — + Np 

to get p - 



3 Diffusion erosion of spectrum 

The process of diffusion leads to destruction of the previous form of spectrum. 
Since the spectrum is rather narrow and (0) takes place one can use initial 
condition as ^-function type. But it is not absolutely necessary. 

There are two approaches which can be used here. At first one can assume 
that the spectrum is narrow according to and use the Green function in 
i/-axis 

G= 7cr xp( — wr- ) 

Here u = pg, 

D = (W + + W-)/2 

W + is absorbtion coefficient and W~ is desorbtion coefficient. These coeffi- 
cients are taken at v = u . 

Another possibility is to take vq as the mean coordinate corresponding to 
z 3 in the previous solution of the relaxation problem. 

The effect of dependence of D on v leads to the following equation 

9 -l = -D u^ 
dt dv 2 
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where 

D = DjvT 

does not essentially depend on v. The last equation can be approximately 
written as 

dj_ &l 
dt °8S 2 

for 

S = z/ 2 / 3 

Then the solution is the Green function and it is known 

N -(S-S„)\ ... 

Here 

C 2 / 3 

Until the half-width 4Dt is many times less than uq there is no difference 
what variable [y or S) is used. 

One can remark that the initial form of the size spectrum will disappear 
rather soon. Really, the characteristic width of spectrum is 

A 2 S = (z + Axrr) 2 - z 2 ~ 2& x xz 

and 

A 2 S _ n A,x 
S ~ Z z 

After the time t 2 the diffusion erosion of the size spectrum will have the 
scale of the initial width of spectrum and the initial form of the size spectrum 
completely disappears. This time allows the estimate 

t 2 = (A 2 S) 2 /(AD) 

Since Aix/ Zfi n 1 one can use both x or S representation of solution. 

One has to note that in the process of relaxation ( varies and D = 
-D(C) varies. The law ((t) is already known from solution of the regular 
evolution equation. Then one can suggest the following approximate account 
of variation of D: instead of Dt in denominator one can approximately use 
JD(((t))dt. 
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One has also to refine the regular evolution equation, which allows to 
prolong it. Really, having written the law of growth 

^£ = ih-K) 

dt t z' 

with 

i 

This differential equation will be more accurate and complex. The problem 
is that (if begin to depend on t and z. 

But here fortunately appears the monodisperse approximation. Since 
Aix <C Zfi n then 

z 3 ^ « Nz 3 > z 2 ^ « NA lX z 2 > z^f « JVAfc* » ^ « JVAfc 

So, one can use the monodisperse approximation. The monodisperse approx- 
imation is rather accurate. In this approximation 

( = $ - Nz 3 

Then 

dz _ $ - Nz 3 2a 

dt ~ t ^ ~ 3<5>z-3Nz 4 ' 
The main momentum is The main advantage is that this momentum 
doesn't change in the real regular law of growth instead of the asymptotic 
one. The last equation can be easily integrated 



/ 



dz t 



(*-^ 3 )(l-3l^) - 



Since < dz/dt < (/r for z > z c one can see the "narrowing" of the 
spectrum. 

Now we have to extract the limits of a "pure diffusion". One has to take 
into account that in reality the pure diffusion equation takes place only in 
the "extended near-critical zone". While the boundary of the "near critical 
zone" is determined by condition 

F c - F ~ 1 
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the "extended near-critical zone" (ENZ) is determined by two conditions 



du ^ 

u<(2 + 3> c = P 3 u c 

The second condition is necessary for the case ( C 1 to exclude the infinite 
zone of big u. If we require that the diffusion length AS ~ y/4Dt has to be 
compared with the regular length AS ~ ( 2 t 2 /T 2 , then we come to p 3 ~ A/(, 
which is less than p c = 2a /3(. 

Certainly, since now dz/dt depends on z the values of pi will change. So, 
here there are only two possibilities: 1). to consider pi slowly varying; 2). to 
solve everything in the monodisperse approximation. 

One can see that both "near-critical zone" and "extended near-critical 
zone" are not symmetrical in respect to u c . This requires to use more so- 
phisticated methods instead of the leading term of the steepens descent to 
calculate the Zeldovitch' factor. But this can be considered as technical 
details. 

The lower boundary where the first condition fails even for small ( de- 
pends on the concrete value of surface tension. We denote the boundary 

as 

ui = u c a pi = p c a 1/3 

The value of a can be estimated when we require that the derivative of 
the free energy equals by the absolute value to the derivative of the free 
energy for big supercritical embryos. This leads to a — 1/8. 

The spread of the gaussian occurs up to the moment 

H = {l-a 2 / 3 fSl/{AD Q ) 

where S c corresponds to the final coordinate. 

When the droplet attains u\ then rather soon (in comparison with the 
time of attaining au c ) it will be dissolved. The probability of growth back is 
small. 

When the essential part of spectrum is in the ENZ one can use the pre- 
sented solution. One can observe and prove the relaxation to the quasi- 
equilibrium. In quasiequilibrium the behavior of p c is regulated by the bal- 
ance condition 

/ f(p,t)3 P 2 f t dp = 
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(f(p,t) is the known distribution) or by 

J r p 

It is important that the vapor consumption by droplets with a known spec- 
trum of sizes regulates the behavior of p c contrary to the method of Lifshitz- 
Slezov where the behavior of p c regulated the form of spectrum. 

Actually nothing happened when f(p,t) is narrow. Let it be like 5- 
function S(p a ). Then p a relaxes to p c and the process stops until the diffusion 
makes spectrum wider. 

4 Further evolution of spectrum in a near 
critical region 

On the base of solution (JHJ) one can see the transformation of spectrum. The 
spectrum is like a Gaussian and according to the balance condition p c is near 
the head of the Gaussian. Then it is clear that at first the lower boundary of 
ENR will be attained by the essential part of spectrum. The upper boundary 
can not be attained first. 

We see that the essential part of spectrum is localized in the region where 
| v — v d <C 1 . Then at first there is no difference whether to use the diffusion 
equation in i/-scale or in S-scale. 

The time when the essential part of spectrum attains p\ can be easily 
calculated 

Wt = {y c - v { ) 2 

Here we calculate time from the moment of the end of the previous period. 

After this moment we have to take into account that practically all 
droplets at pi are going to dissolve. Then 

= o 

Solution of diffusion equation (it has to be solved now in S'-scale) with 
such a boundary condition is 

f(p) = o 

for 

P < Pi 
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and 

f(S) = G + - G_ (4) 

for 

S 1/2 > Pl 

Here 

N (S-S )\ 

(s+ — i ex P( 7T; ) 

N (S-S 1 )\ 

and So and Si are defined by conditions that So = z 2 in and S — Si = Si — Si 
(where Si = v 2 ^). 

A special question appears what is here D. But at first we shall analyze 
the properties of solution. One can see that contrary to the unbounded 
solution with a half-width 

AS ~ t 1 ' 2 

and the maximum (and mean) value is staying 

Srnax = COHSt 

the new solution has the moving maximum value 

Smax ~ t ^ 

the moving mean value 

Smean ~ t ^ 

and the growing half-width 

AS < t 1 ' 2 

which ensures the relative localization of spectrum. 
So, we have the following estimate 

A < <7 

The value of supersaturation has to be determined from 

/ o °°^3p 2 /(p)f = (5) 
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Then we get p c and the supersaturation. 

Then in this solution we have to take into account that D « D(p c (Q). 
Again we can approximately substitute Dt by / D(((t))dt. 

Actually, this is the end of solution. Everything presented later refers to 
rather small tails which can be seen only in enormously big systems. 

5 Formation of a tail 

The previous equation to get the supersaturation is not a unique way to 
do it. Rigorously speaking the supersaturation has to be derived from the 
substance balance equation 



Distribution has at first to be taken here from 

One can see that the part of the size spectrum with big p is important 
in the substance balance due to the factor p 3 . But the law of growth of 
the droplets of big sizes is regular, it is not the diffusion walking on the flat 
potential surface as it was supposed in derivation of solution 

We suppose that at 



the law of growth will be changed to the regular asymptotic law of growth 



dt t 

Here (3 is a parameter, the normal value is 

(3 = 1.5 

The rate of appearance of droplets at the beginning of the region of regular 
growth can be calculated as following. At first we can calculate f(p) by 





Pr = f3p c 



dp C 



f(p) = f(S) 



dS 



f(S)2S 1/2 



dp 



on the base of solution (J3J). Then 



J = f(M 



C 



T 
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or 

J= f(P Pc )-(l-P) 

T 

for the intensity of formation of the tail. Then we get the tail formed accord- 
ing to 

ftail(p,t) = f(P P c(t'))^ 

where t' is defined as 

p = (3p c + f 



T 

Let z be a coordinate p of some droplet 2 . Then x is z — p for the given 
droplet. It remains constant during the regular growth. Then 

ftail\Pi t) = ftail\ x ) 

and it is a known function. 

One can see the following approximation. Certainly, the position of S c 
coincides in initial moments of time with So and do not vary essentially. It 
a quasi-integral of evolution. Then denoting 

2a 

Co — - — 

3p 

one can put Co instead of C(t") in the previous formula. The movement of the 
boundary S r = (/3p c ) 2 can be neglected and one can act in a quasi- stationary 
approximation. As the result one gets the form of the tail. 

More accurate is to take S c from the solution (j3J) ignoring the formation 
of the tail. Then 

C ° (t,,) = WW) 
Now we can present the distribution as 

/ = f± + ftail 

and in the quasiequilibrium we have 

J o Ca dp3p 2 (f ± + f taU )^ = (7) 

This gives p c in quasiequilibrium and a new value of supersaturation. This 
closes the iteration loop. Actually only one loop is necessary. 

2 Morc convenient is to choose z as the coordinate of the front side of the tail. 



11 



6 Dissolution of the head of spectrum 

At first the tail is not important in the substance balance. But the quantity 
of substance G ta u in the tail grows faster than something proportional to t 3 . 
So the droplets in the tail become the essential consumers of vapor. 
Now the distribution density is the following 

f(p) = o 

for 

P< Pi 

and 

f(S) = G + - G_ (8) 

for 

f ftail 

for 

P> Pr 

Then one can get the supersaturation from the balance equation 
C = <£>* - / f(S)S^ 2 dS + / f tml (p)p 3 dp 

J Si J p r 

At first the vapor consumption of the droplets in the tail i.e. the evolution 

G t ail(t)=J s ftail(P, t)dp 

determines the value of supersaturation 

de~ dt tail 

The head of the spectrum begins to dissolve since the droplets in the head 
of spectrum become mainly subcritical. This diminishes the quantity 

Ghead = / f(P,t)dt 

and partially stabilizes the supersaturation. 
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The result will be the dissolution of the head of the size spectrum. The 
final of this period will be when 

Ntail ~ Ntotal — N tai i 

One can propose the approximate (not necessary) model: at first the 
intensity of vapor consumption by the tail can be calculated according to 

i Gtau = 3 C ftMp2 tt dp 

Then the number N head in the head satisfies the following equation 

[jN head \4 + jGtail = 

The growth of G ta u occurs in the avalanche manner. So, the dissolution of 
the head occurs in the avalanche manner also. So, this solution can prolonged 
up to 

Ntail < N - N tail 



7 Dissolution of the essential tail of spectrum 

One can spread the solution to the beginning of the dissolution of the tail of 
spectrum. 

But when the complete dissolution of the head is finished one can propose 
the more effective procedure. 
From 

where the changing momentums are 

/X 
X l ftail(x)dx 
-oo 

(here x corresponds to p = /3p c ). are known functions of x c we take x c as 
a function of z. The region p < f3p c is negligible in definition of pi when 
z » p c . Now the law of growth 

dz 2a 
dt 3(z — x c ) 
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becomes the closed differential equation of the first order. Since the r.h.s. 
of the last equation is known function of z the last equation can be easily 
integrated 

f dz 
t = r J —is— 

3(z-x c (z)) 

It can be rewritten as 

i 

Since /ii(x(z)) can be well approximated as polynomials on z, the integration 
can not cause difficulties. 

Here one can see the important property which allows to solve this equa- 
tion very simply: the vapor is consumed by the droplets with regular growth. 

One can show the simplicity of solution by the following example: suppose 
that the height of tail /* is constant. Then 

G ta a = f*(z 4 - 43 4 ) 

and 

d7 

Tf t =$-Mz 4 -zj(z)P*) 

One can integrate this equation. Rewrite this equation for evolution of z c (t) 

^ = $ + U*\P - 
3z c 

Then z can be expressed through z c 



z = /,- 1/4 ($ + u£P - f^) 174 



Having differentiated the last relation and substituted it into the initial bal- 
ance equation 

*-/.(* 4 -/.*^)=T- 



1 dz 4 



3z 3 dt 
we get 
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which can be integrated 

t = rj ^f; l '\fMP 4 + + - ^-r m dz c 

The end of this period takes place when the number of droplets in ENR, 
namely N ENR becomes comparable with the rest droplets 

Nenr ~ N tota i — N ENR 



8 Closure of the loop of evolution 

Later one has to take into account the diffusion of the rest of the tail. Again 
we have the spectrum of some form in ENR. Since the form of the tail is 
rather flat we take the flat spectrum with a gaussian front side. At first 
we see the dissolution of the flat part of the size spectrum. It is easy to 
describe because we see the one dimensional heat conductivity problem with 
a boundary condition of a following type: at S — Si the distribution P is 
zero; at big S the distribution P vanishes, then at S = z 2 one can see the 
free diffusion. The initial condition is 

P( x , t = 0) = Q(S - Si)AQ(S k - S) 

where A is the amplitude of tail (constant or slowly varying function) and 
S k is the position of the front side of tail (in shifted coordinates), 6(x) is the 
Heavisaid's function. 

To keep condition P — at S — Si we can prolong the previous initial 
condition oddly. 

The solution of diffusion problem can be found anywhere. 

Having moved in S'-scale the origin to the former Si we see the solution 

P = J F(x,0)exp(- { -^f)dx 

Here F(x, 0) is the distribution given by 

F(x, t = 0) = Q(S)G(S k -Si- S)A - Q{-S)Q{S + S k - S{)A 

It is quite clear that when A is constant or can be approximated by poly- 
nomials then the solution can be written via a linear combination of error 
functions. 
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It is also clear that when 

t » Sl/{AD) 

then the solution begin to resemble the gaussian. 

When the flat region is dissolved the size spectrum resembles the gaussian 
and we come to the already studied problem. So the procedure of all sections 
starting from the diffusion erosion of the size spectrum in the ENR has to be 
repeated. This closes the loop and these loops will be inevitably continued 
in future. One has to stress that the number of droplets radically decreases 
in every loop and to see these loops one has to observe an extremely giant 
system. To give description on the base of continuous approximation even in 
the beginning of the tail of the size spectrum there have to be many droplets. 
Earlier or later diffusion and stochastic growth of droplets lead to the change 
of the regime of evolution presented here. Also the change of regime of vapor 
consumption, the heat release effects, coagulation etc. destroy the presented 
solution. 
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